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Abstract 

We consider random wallcs on the surface of the sphere Sn~i {n > 2) of the n-dimensional 
Euclidean space En, in short a hypersphere. By solving the diffusion equation in Sn-i we 
show that the usual law (r^) oc t valid in E^-i should be replaced in Sn-i by the generic law 
{cos9) oc exp(— t/r), where 9 denotes the angular displacement of the walker. More generally 
one has (C^^^ ^ cos(^)) oc exp(— t/r(L, n)) where C^^"^ ^ a Gegenbauer polynomial. Conjectures 
concerning random walks on a fractal inscribed in Sn-i are given tentatively. 

PACS numbers: 05.40.Fb,05.40.Jc 
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I. INTRODUCTION 



We consider a random walk (RW) on the surface of a sphere Sn-i of the real Euclidean 
space En of dimension n, a hypersphere for short. Let O be the center of the hypersphere 
and R its radius, Sn-i is thus defined as the set of points M = (xi, . . . , x„) of En such that 
^2^=1^1 = R^- If the legs of the walker are much smaller than the radius R of the sphere 
then the diffusion process can be described by a continuous diffusion equation, i.e. 



where p{M,t) is the density probabibily for the walker to be at point M at time t. D{> 0) 
is the diffusion coefficient and A'^""^ denotes the Laplace-Beltrami operator in space Sn-i- 
Explicit initial conditions will be specified later. Note that for imaginary times eq. (jl.ip 
identifies with the Schrodinger equation of a free particle of Sn-i-, alternatively this equation 
can be considered as describing the diffusive rotational motion of a n-dimensional linear 
polar molecule. 

The solution of eq. (jl.ll) in the case n = 3 is known since the work of Debye on Brownian 
rotors.i^ In this paper we extend Debye's proof to arbitrary dimensions n and show that a 
salient feature of a RW on a hypersphere is the validity of the generic law 



where the unit vector ^ = OM/R denotes the orientation of vector OM in En- With the 
north pole of Sn-i chosen to be M(t = 0) the angle 9 in eq. (jl.2|) is thus the colatitude of point 
M in spherical coordinates. The brackets (. . .) denote a spatial average with weight p(M, t), 
and the dot in the LHS is the usual scalar product in En- The result fll.2|) . given without 
proof in ref. O, is not surprising; indeed as noted by Jund et al. when Dt/R^ <^ 1 one can 
expand both sides of eq. ()1.2j) yielding for the projection x of vector OM in the euclidean 
plane En-i tangent to Sn-i at point Mq to the behavior {x'^) ~ R^{9'^) ~ 2t{n — 1)D 
characteristic of a random walk in an euclidean space of {n — 1) dimensions. 

Recently Nissfolk et al. have given a solution of p. HI in the special case n = 4.— Their 
expression for p(M, t) is different from that which we derive in sec. IIIII and, in particular, 
does not allow to recover easily eq. ()1.2|) . The equivalence of the two solutions is established 
in the last section of this paper. 




(1.1) 



(^(0) ■ m ^ (COS0 (t)) = exp {-Din - l)t/R') , 



(1.2) 



Our paper is organized as follows. In next section |n] we give an elementary proof of 
eq. ()1.2j) . The expression of the Green's function of eq. (jl.lj) is derived in section UTTl from 
which eq. (jl.2|) can deduced. The solution, inspired from the seminal paper of Debye on 
3D rotors^'^, requires the whole machinery of hyperspherical harmonics and the use of some 
properties of the rotation group SO{n)~. Finally in section HVI we consider in more details 
the case n = 4. Our solution of for n = 4 and that proposed by Nissfolk et alA are 
shown to be identical. A preliminary account of the present work, however devoted to the 
special case n = 4, can be found in ref. Q. 

II. AN ELEMENTARY PROOF OF {cosO) oc exp(-t/T) 

We suppose that the walker it standing at some point OMq = R^q of Sn-i at time t = 0. 
At subsequent times its dynamics is governed by eq. (jl.ip and we study the time behavior 
of the random variable cos 6 = ■ ^- We choose the north pole of Sn-i to be Mq and 6 is 
thus the colatitude of point M in spherical coordinates. We define 



{cose){t)= p{M,t) cose dT{M) (2.1) 

where dT{M) is the infinitesimal volume element in Sn-i- It follows from eq. (jl.ll) that 

d 



(2.2) 



(cos^) = / D cose A''"-'p{M,t) dT{M) 
= f Dp{M,t) A^-' cose dT{M) 



where we have made use of Green's theorem in Sn-i^ At this point we recall that the 
Laplacian A^" in En may be decomposed as the sun>^ 

A-^" = A^ + A-^"-! , 

where A^ acts only on the variable R and A"^"~% which acts only on angular variables, is 
called the angular part of the Laplacian or, alternatively, the Laplace-Beltrami operator.— 
From these remarks, the action of A"^"-^ on cos^ will be easily obtained. In one hand we 
have obviously A^"{OM ■ ^o) = 0; however, on the other hand, it follows from eqs. ()2.3p 
that A^"{OM ■ ^o) = A^"{Rcose) = (n - 1) cose/R + RA''^"-^ cose whence A"^"-! cos 6* = 



— {n — 1) cos9/R^ yielding, when inserted in eq. ()2.2|) . to the simple equation 

^(cos^) + {n- l)^(cos^^) = , (2.4) 
the solution of which being of course given by eq. ()1.2|) . 

III. THE GENERAL SOLUTION 

We shall make extensive use of the properties of the hyperspherical harmonics to solve 
eq. (Il.lj) . A short resume on these functions might prove useful for the reader. We will 
adopt the definitions and notations of the chapter IX of the classical textbook by Vilenkin.— 
In the case n = 4 this yields to a definition of the harmonics which is slightly different from 
that which we used previously.^*^*^ 

To avoid confusion we shall denote by the hj^ersphere of radius R = 1 (When, 
incidentally, we talk about the volumes of the spaces S'„_i or iS„_i we will mean in fact the 
areas of these spheres considered as manifolds of En). The spherical coordinates of the unit 
vector ^ = OM = (^i, . . . , ^„) of Sn-i are defined by the relations 

^1 = sin 6'„_i . . . sin 6^2 sin 6*1 , 

,^2 = sin 6n^i ... sin 612 cos 6^1 , ... , 

= sin6'„_i cos6'„_2 , 

= cos6'„_i , (3.1) 

where < < 27r and < 6k < it for k ^ 1. Note that 6n~i is the colatitude of point M 
and was denoted by 6 in section |TT1 The integration measure in iS„_i will be defined as 

d$, = sin"-2 9n-i ... sin 62 dOi . . . dOn-i , (3.2) 

where An-i = 27r"/^/r(ra/2) is the surface of the sphere Sn-i- With this normalization 

d^ = l , (3.3) 

<Sn-l 

and the infinitesimal volume element in Sn-i is given by dr = The expression 



of the Laplacian A*^"-! = ^Sn-\ spherical coordinates reads as 

sin^ On-i sin""^ 6'„_2 <96'„_2 ""^ <96'„^2 



sin' e„_i...sin2 62 dOf 

The (hyper) spherical harmonics '^l,k{^) are defined to be the eigenvectors of the operator 
^<s„_i_ rjj-^g spectrum of A'^"~'^ is given by the integers X{n,L) = —L{L + n — 2) where 
L = 0, 1, ... is a positive integer. One can show that K = {ki, . . . , ±A;„_2) with L > ki > 
^2 > kn-2 > 0. Therefore there are exactly h{n, L) = (2L + n — 2){n + L — 3)1/ {n — 2)\L\ 
distinct harmonics '^l,k{^) corresponding to the eigenvalue X{n,L). The H2,,j<:(^) constitue 
a complete basis set for expanding any square integrable function / G C^[Sn~i] defined 
on Sn~i- An explicit expression of "Bi^Ki^) will be of little use and can be found in the 
Vilenkin.— As the usual spherical harmonics the H2,,k:(^) satisfy to the following properties*^ 

• (i) orthogonality 

f EIj,{^) El',k'{^) = 6l,l' 6k,k' , (3.5) 

Sn-l 

• (ii) completeness 

E^U(0 2l,k(O = 5^"-H^,O, (3.6) 

L,K 

where 6'^"~^ is the Dirac distribution for the unit hypersphere defined as 

/ /(0 5''"-H^,Orf^ = /(O- (3.7) 

• (iii) addition theorem 

E-i-(^) ^^Me) = ~ ' cr~\^ ■ a , (3.8) 

K 

where the dot in the RHS denotes the usual scalar product in En, i.e ^ ■ ^' = cosip 
where ip is the angle between the two unit vectors $, and and C^^^"^ is a Gegen- 
bauer polynomial. The Gegenbauer polynomials are a generalization of Legendre 
polynomials. is defined as the coefficient of in the power-series expansion of the 
function 

00 

(l-2t/i + /i2)-?' = ^C];(t)/i^ . (3.9) 

L=0 



We have now in hand all the tools to solve eq. (Il.lj) . We first expand p{M,T) in terms of 
the El,k, 

oo 

p{M,t) = Y,Y.P^'^^^) -^.^(^) (3.10) 

L=0 K 

and then insert the expansion ()3.10|) in eq. (jl.lj) . Making use of the orthogonal properties 
of the Sl^x yields an infinite system of non-coupled equations 

d L(L + n-2)\ , , 

' D pLMt) = , (3.11) 



dt i?2 
the solution of which reads obviously as 

PlM^) = Pl,x(0) exp {-DL {L + n - 2) t/R^) . (3.12) 

We shall denote by p{M, t|Mo, 0) the solution of ()1.1|) corresponding to the initial condition 
p(M,0) = 5^-1 (M, Mo) = 5'5"-i(^,^o)/^""^A„_i, i.e. the solution of 

Pp(M,t|Mo,0) = 6it) ^^J^^ ^ . (3.13) 

It follows readily from eqs. ()3.12|) and ()3.6|) that the Green function p{M,t\Mo,0) can be 
expressed as 

p(M,t|Mo,0) = (t < 0) 

^ oo 

p{M,t\Mo,0) = E E -U(^") -^^^(^) 

^ ^^-^ L=Q K 

X exp{-DL{L + n-2)t/R'^) (t > 0) . (3.14) 



Eq. ()3.14|) can be further simplified with the help of the addition theorem ()3.8|) yielding our 
final result 



" ^ L=0 

X exp{-DL{L + n-2)t/R^) (t > 0) . (3.15) 



Some comments on eq. ()3.15j) are in order. 

(i) The solution ()3.15|) is invariant under rotation about the axis a^s expected. Eq. 
()3.15|) is a generalization for all n of Debye's result which corresponds to the case n = 3.^^ 
In this case the Gegenbauer polynomials reduce to the Legendre polynomials P^. In the 



case n = 4 we recover the result of ref. Q where the Gegenbauer polynomials reduce to 
Tchebycheff polynomials of second kind. To be more precise, recall thalj^ 

Cl^\cose) = Piicose) {n = 3) (3.16) 
Cl(cos«) = 5Hi^ („ = 4) (3.17) 

(ii) Note that ^ p{M, t|Mo, 0)dT = po,o = 1, i-e. the probability is conserved, and that 
we also have limt_++oo p{M, t|Mo, 0) = 1/Qn-i, where Qn~i = -R"~^^n-i is the volume of the 
space Sn— 1, i.e. the solution of the diffusive process is uniform after infinite time. 

(iii) Let us define the time correlation functions 

= [ di, I di Cl'^-\i,-i)p{M,t\M,,Q), (3.18) 

1 J Sn—\ 

where we have introduced the reduced Green's function p = i?"~^A„_ip. As a consequence 
of the orthogonality properties ()3.5j) we find that 

:Fr-\t) = ^I'^^^f exp {-DLiL + n- 2) t/R^) (3.19) 

Special cases are of interest. Firstly, since for L = 1 we have Cf"^^ ^(-u) = (n — 2)u, equation 
(j3.19|) indeed gives back eq. (jl.2|) . For n = 3 and 4 we recover the results for the 3D and 
AD rotorsii^ii 

{Pl (cos 6)) = exp {-DL {L + 1) t/R^) {n = 3) (3.20) 

smCL + l^^ = (L + 1) exp (-DL {L + 2) t/R^) (n = 4) , (3.21) 
sm 6 

Defining now the reorientational time tI^"^"^ as 

r?/^ ^ = / / ^ = , 3.22 

we have now at our disposal a Kubo formula for the diffusion coefficient D. Note the 
aesthetic relation t^"^'^ jr^i'^'^ = L'{L' + n — 2)/L{L + n — 2) which generalizes Debye's 
result to arbitrary dimensions.-^ Since simulations of real 3D liquids or plasmas are feasible 
(moreover efficient) in S^^'^^ eq. ()3.22|1 is of prime importance since it should allow the 
computation of the self-diffusion coefficient of such systems in the course of equilibrium 
molecular dynamics simulations. 



(iv) Simulations of random walks on fractals inscribed in Sn-i were reported recently.— 
While it is known that for random walks on fractal clusters in Euclidean spaces the anoma- 
lous diffusion law cx hold o^^i^^ (where f3 is some ejcponent depending of the 
fractal dimensions of the cluster) the numerical results of ref. |3| give evidence of the law 
{cos6) oc exp(— (t/r)^), i.e. a stretched exponential relaxation, on the hypersphere. As sug- 
gested by Jund et al^ the anomalous diffusion law for a RW on a fractal can be understood 
by replacing the time t by a fractal time in the diffusion equation and thus in the 
solution ()3.15|) . With this assumption we now have the ansatz 

[rT^'/rT^') = L'{L' + n- 2)/L{L + n - 2) , (3.23) 
which could be checked in numerical simulations. 



IV. RANDOM WALKS IN S-^ 

We specialize to the case tt, = 3 and want to show that the expression ()3.15|) of 
p(M, t|Mo,0) is equivalent to that of ref. jj]. This can be done as follows. Let us rewrite 
the reduced p = 27:"^ p{M,t\MQ,Q) as 

sin(L + l)^ (4.1) 
^ — ^ sm-w 

L=0 ^ 

where K = Dt/B?. A priori the angle if) is in the range (0, vr), however since the function 
is formally even in ip we define ^{—ip) = 'p{'ip) for negative angles. This gives us a periodic 
function of period 2tt defined for all ip We introduce now the periodic function 

ri' 

F(V') = / d^' p{tp') sin(V'') (4.2) 







which can be rewritten after some algebra as 



X K 

F(^) = Fo - Yl i-^P^) exp(-^pV') (4.3) 



2 

p=—oo 



where Fq is some unessential constant independent of angle ip. At this point we recall 
Poisson summation theorem which states that for any function ip{x) holomorphic in the 
strip —a < < a one has 



+ 00 +CO „ 



^{y)e'Py dy . (4.4) 



p=— oo 



Applying Poisson theorem for the Gaussian we get 



FM = - J: exp [-'-^j^ ) . (4.5) 

n=— oo ^ ' 

which, after differentiation yields for p 



which coincides with the result of ref. y, apart the prefactor which is not specified. Relations 
similar to eq. ()4.6|) can be obtained for even values of n (other than 4) but we failed to get 
anything similar for odd n's. 
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